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The pressure dependence of the inverse square of the magnetic penetration depth A~ 2 
in k-(ET) 2 Cu(NCS)2 was measured by Larkin, et al. According to the paper, A -2 behaves 
differently under low pressure and under high pressure. Under low pressure, the development 
of A~ 2 just below T = T c is rapid compared to the case under high pressure. Moreover, T c in 
k-(ET) 2 Cu(NCS)2 increases under c-axis pressure up to lkbar and decreases under higher 
pressure, while T c decreases monotonically under the hydrostatic pressure, or under the 
uniaxial pressure parallel to other axes. In order to explain these behaviors, we calculate T c 
and A -2 for k-(ET)2Cu(NCS)2 under pressure. In the calculation we mainly use an effective 
dimer Hubbard model. In conclusion, the behavior of A~ 2 results from three effects: the 
variation of the bandwidth of quasiparticles, the change of the Fermi surfaces, and the effect 
of vertex correction. This is a different mechanism from that of A~ 2 in cuprates which we 
observe when the doping varies. Moreover, we explain the increase in T c under the c-axis 
pressure up to lkbar and the decrease in T c over lkbar from our calculation. With the 
increase in the c-axis pressure, two competitive effects with respect to T c appear. One is 
the approach of the Fermi surface to the antiferromagnetic Brillouin zone boundary, and 
the other is the suppression of the electron correlation. Under the low c-axis pressure, T c 
increases since the former effect is dominant. On the other hand, T c decreases since the latter 
effect is dominant under the high c-axis pressure. 

KEYWORDS: unconventional superconductivity, organic superconductor, penetration depth, 
vertex correction, pressure 



1. Introduction 

Recently, many measurements have been performed in the superconducting state in or- 
ganic conductors including /i-(ET)2Cu(NCS)2- Many interesting behaviors are observed under 
pressure in organic conductors, since they are soft materials. Among them, we study the pres- 
sure dependence of the inverse square of the magnetic penetration depth A -2 and the behavior 
of the superconducting transition temperature T c under the hydrostatic and the uniaxial pres- 
sure by using the effective dimer Hubbard model. 

We calculate A -2 in k-(ET)2Cu(NCS)2 on the basis of Jujo's theory. 1,2 This theory is 
a general theory on the transport phenomena and the magnetic penetration depth in the 
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superconducting state, including the effect of the vertex correction, and can be applied to 
unconventional superconductors such as cuprates and organic superconductors k-(ET)2X. 
By taking account of the vertex correction, which can be considered as the backflow effect, 
this theory explains why the value of inverse square of the magnetic penetration depth A -2 
is suppressed near T = under the strong antiferromagnetic fluctuation in unconventional 
superconductors, such as cuprates in the underdoped region. In this paper we discuss the 
pressure dependence of the extrapolated value of A~ 2 at T = and the development of 
A~ 2 just below T c . Moreover, T c under pressure is calculated by the linearized Dyson-Gor'kov 
equations, and the different dependence of T c on hydrostatic and uniaxial pressure is explained. 

2. The Pressure Dependence of A -2 and of T c in k-(ET) 2 Cu(NCS) 2 



The pressure dependence of the inverse square of the magnetic penetration depth A~ 2 in 
K-(d 8 ET) 2 Cu(NCS)2 was measured by Larkin, et al. 3 Here, we refer to k-(ET) 2 Cu(NCS) 2 in 
which hydrogens in ET molecules are deuterated as K-(<i8ET)2Cu(NCS)2, while the materials 
with protonated ones are referred to as K-(/igET)2Cu(NCS)2. Ref.3 shows that the value of 
is large under low pressure (P = Obar) compared to that under high pressure(P is 
up to 1290 bar). This implies that there is a rapid development of A -2 just below T c under low 
pressure. Moreover, Ref.3 shows that A~ 2 is suppressed in the low temperature region under 
low pressure, although the data near T = is insufficient. This suppression is analogous to 
the doping dependence of A~ 2 in high-T c cuprates such as LSCO. 4 The A" 2 -T diagram in the 
underdoped and overdoped LSCO behaves similarly to k-(ET)2Cu(NCS)2 under low pressure 
and under high pressure, respectively. The behavior of A -2 in cuprates is understood from the 
fact that the effect of the vertex correction is strong in the underdoped region, and weak in 
the overdoped region. The vertex correction affects strongly under low pressure and weakly 
under high pressure in /-c-(ET)2Cu(NCS)2. In this paper, we study the behavior of A~ 2 in 
k-(ET)2Cu(NCS)2 theoretically by carrying out the actual calculation. 

Next, we refer to the pressure dependence of T c for k-(ET)2Cu(NCS)2- T c under the 
hydrostatic pressure have been measured by many experiments. 5-7 According to the results, 
T c decreases monotonically as the pressure increases. On the other hand, the behavior of 
T c under the uniaxial pressure is quite different from that under the hydrostatic pressure 
according to the recent measurements for k-(ET)2Cu(NCS)2. 8 ~ 10 When the uniaxial pressure 
is applied in parallel to 6-axis or to a-axis, T c decreases. This is similar to the behavior 
under the hydrostatic pressure. On the other hand, T c increases up to lkbar under the c-axis 
pressure and decreases under higher pressure, which is in contrast to the behavior of T c under 
the hydrostatic pressure. By calculation, we show that this is because two competitive effects 
with respect to T c appear under the c-axis pressure: the approach of the Fermi surface to the 
antiferromagnetic Brillouin zone boundary and the suppression of the electron correlation. 
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Fig. 1. (a): The original lattice structure of /c-(ET)2X, where the ellipses represent ET molecules. A 
dimcr consists of two parallel ET molecules, (b): The lattice structure of k-(ET) 2 X when the pair 
of ET molecules are dimerized. 11, 12 The circles represent dimers, and the rectangle represents 
the unit cell. k-(ET) 2 X is quasi-two-dimensional (Q2D) conductor where the layer is usually 
denoted as the &c-plane. x- and y-axes are introduced for the convenience of calculation, a-axis is 
perpendicular to the frc-planc. 



3. Effective Dimer Hubbard model for Organic Superconductor k-(ET)2X 

We use the effective dimer Hubbard model for k-(ET)2X. 5 ' 6 ' 11,12 In this paper we adopt 
a quasi-two-dimensional (Q2D) band structure ignoring the interlayer hopping term for k- 
(ET)2X and symmetrize two dimers in the unit cell. 6 Then the Hamiltonian is expressed as 
follows: 13 

k k,k',q 

= — It (cos k x + cos ky) — 2t' cos(k x + k y ), (2) 

where t and t' are the nearest-neighbor and next-nearest-neighbor hopping terms shown in 
Fig.l. U is the Coulomb repulsion and is the dispersion. The bandwidth W is expressed as 
the maximum of the dispersion measured from the minimum: Eh max — Ekmin ~ 8i. The exact 
value of W is 8.22t and 8.39t at t'/t = 0.7 and 0.75, respectively. In k-(ET) 2 X, the electrons 
are half-filled. In this paper, we adopt the FLEX approximation, the self-consistent second 
order perturbation theory (SC-SOPT), and the third order perturbation theory (TOPT) in 
the calculation of the self-energy terms. 1 ' 2 We adopt different approximations in the calcula- 
tion for A -2 and for T c . In the calculation for A~ 2 FLEX and SC-SOPT are adopted, since 
conservation approximations are needed. FLEX is appropriate for the case with the strong 
antiferromagnetic fluctuation, while SC-SOPT is used to show the weak correlation effect 
and the weak antiferromagnetic fluctuation qualitatively. On the other hand, SC-SOPT is not 
enough to give actual T c since the third order term in U is important for calculating T c . There- 
fore when the antiferromagnetic fluctuation is weak, we adopt TOPT instead of SC-SOPT in 
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the calculation for T c . As a conclusion, we adopt FLEX and TOPT in the calculation for T c . 
The former is used for the strong antiferromagnetic fluctuation and the latter is used for the 
weak correlation. 

Within FLEX and SC-SOPT, the self-energy term and the Green's function in the normal 
state S n is expressed as follows: 

^n(k) = ^^V(q)G(k-q), ( 3 ) 

q 

G ^ = i£ B -( efc - M )-Sn(fc)' W 

where k = (k, e n ), e n = (2n + 1) ttT is the fermion Matsubara frequency, and {i is the chemical 
potential. The normal vertex V (q) in Eq.3 is given by 

for FLEX and 

V(q)=U 2 X (q) (6) 



for SC-SOPT. The spin susceptibility \ (<?) is expressed as 

T 

N 



q 

Note that q = (q,uJi), and uji = 2nirT is the boson Matsubara frequency. Within FLEX and 
SC-SOPT, we obtain the normal self-energy term and the normal Green's function by solving 
Eqs.3 and 4 self-consistently. 

Within TOPT, the normal self-energy term is expressed as follows. 

S„ (k) = ^ i u2 Xo ( k ~ k ') + U *Xo ( k ~ k ') + (k + k')] G (k'), (8) 



N 

V 



N 

k 

k 

where xo (q) an d 4>o (q) are calculated from the bare Green's function Go (k) = ig _^ £k _^ and 
Xo (q) x(q)- This shows that the calculation of the self-energy term and the Green function 
within TOPT is not self-consistent in contrast to FLEX and SC-SOPT. This is because self- 
consistent TOPT seems to be complicated in actual calculation. 

Here we introduce the FLEX approximation in the superconducting state, the normal and 
anomalous self-energy terms S n (k) and E a (k) are given by 

E n (fc) = ^XV»(<z) G (*-?)' (11) 
q 
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where V n (q) and V a (q) are the normal and anomalous vertices in the superconducting state, 
respectively. The normal and anomalous Green's functions G(k — q) and F (k — q) are written 

as 

ie« + (e-fc - n) + S n (-k) 



G{k) 
F(k) 



[\e n - (e fc - n) - S n (fc)] [ie n + (e fc - //) + E n (-AO] - S a (A;) 2 ' 

-S a (k) 



[is n - (e fc - //) - S n (fe)] [ie n + (e fc - //) + S n (-AO] - S a 
V n (q) and K (g) are given by 



U 2 



x s (q) 



Xs {q) + 



Xc{q) 



V a {q) 



2 L l-Uxs(q) 
U 2 \ Xs (q) . , 



1 + U X c (q) 

Xc(q) 

1 + U Xc (q) 



- Xc (q) 
+ Xc (q) 



with 



Xs(q) = -^^ G ( k + ^ G ^ ±F ( k + ^ F ^ 



(13) 
(14) 

(15) 
(16) 

(17) 



where + and — correspond to Xs and Xc: respectively. 

We solve Eqs. 11-17 self-consistently to calculate G (k), F (k), S n (k) and T, a (k). The sym- 
metry of Cooper pair in k-(ET)2X is d x 2_ y 2, with which the superconducting gap E a (k) has 
nodes in ±7r/4 directions in the /c-space. 13 The d x i_ y i gap symmetry for k-(ET)2X is justified 
later in calculating T c . In order to realize this symmetry, we set E a (k, e±i) oc (cos k x — cos k y ) 
as the intial value for solving Eqs. 11-17. 

4. Calculation of T c by Dyson-Gor'kov Equations 

In this section, we calculate T c in k-(ET)2X by the linearized Dyson-Gor'kov equation 
within FLEX 14-18 and TOPT. 13 First of all, the normal and anomalous Green's functions 
satisfy Dyson-Gor'kov equations. 19 

G (k) = G (k) + G (k) S n (k) G (k) + G (k) S a (k) F* (k) , (18) 

Ft (k) = G (-£:) S n (-k) Ft (k) + G (-fc) S a (-k) G (k) . (19) 

When T approaches T c within the superconducting state, Eqs. 18 and 19 are linearized as 
follows, since F (At) < G (k). 

F{k) = |G(A:)| 2 S a (A;), 



G{k) = G (A:) + G (k) S n (k)G(k). 



(20) 
(21) 



Within FLEX, S n (A;) and E a (k) in the vicinity of T c are calculated by linearizing Eqs. 11 and 
12, respectively. Then, S n (A;) is the same as that in the normal state, which is expressed by 
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Eqs.3 and 5, while E a (k) is expressed by the linearized Dyson-Gor'kov equation: 

s - (*) = - S E y - («) i G ( fc - 9)i 2 ^ ( k - 1) > 



k (?) = ^ 2 



3 x(<?) 1 x(<?) 



+ £/. 



L2l-x(?) 2l + x (g). 

On the other hand, the linearized Dyson-Gor'kov equation for TOPT is 



SaW = --y£ V a (k, k') \G (k') | 2 S a (A/) , 



where, 



K (fc, A/) = Vrpa (fc, fc') + V^t {k, k') , 

Vrpa (k, k')=U + U 2 X0 (k + k') + 2U\l {k + *0 , 



(22) 
(23) 

(24) 

(25) 
(26) 

(27) 



Vvert (k, k') = 2U 3 -Re G (h) G (k + h - k') [ X o (A; + h) - 4> (k + h)]. 

ki 

If we replace the left hand side of Eqs.22 and 24 by aS a (k), this equation can be considered 
as an eigenvalue equation with eigenvalue a and eigenvector S a (k). T c is the temperature at 
which the maximum eigenvalue a max reaches to unity. S a (At) with which the largest eigenvalue 
is obtained represents the superconducting gap symmetry. Among several gap functions, the 
d x 2_ y 2 state possesses the maximum eigenvalue in the region 0.17 < t/U < 0.35 and 0.4 < 
t'/t < 0.8 within FLEX, and 0.14 < t/U < 0.25 and 0.4 < t'/t < 0.8 within TOPT. 

5. Effect of the vertex correction and A -2 : Introduction to Jujo's theory. 

In this section, we introduce Jujo's theory shortly. From Kubo formula, the penetration 
depth is expressed with the electromagnetic response kernel as follows. 



A" 2 = -4ttK^ (q - 0) , 



(28) 



(?-0) = -^£>[AJ ^ k + q) G(k + q)K (k + q,k)G(k) 
k,n 



E 



q^o V f-^ dk v dka 

k,n ^ 



Tr 



f 3 G(k) 



j£n0+ 



(29) 



Equation 29 is expressed by Nambu matrices. The Green's function matrix is defined as: 

' G(k) F(k) \ 
F{k) _ G( _ fc) ) 

1 



G(k) 



ie n f - e k f 3 - £ (k) 
where the self-energy matrix £ (A;) is given by 

' E n (k) £ a (fc) \ 
v E a (fc) -E n (-fc) J 



(30) 



t(k) 



(31) 
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Matrices fo and T3 are given by 

-(::)• 

A° (*,* + «) =( Vk ; ° ) (34) 



is the bare three-point vertex, where = ^ is the bare velocity of electrons. A M (k + 5, k) 
satisfy the following integral equation: 

+ ^y2y2 r ^ j (k + q,k'-k' + q,k)(G(k' + q )A^(k' + q ,k / )G(k')) , (35) 

iv — — V / m,n 

A;' m,n 

where irreducible four-point vertex Ti n ^ m j is expressed by the functional derivative of the 
self-energy matrix by the Green's function matrix: 

G 



F . - 



(36) 



After rather long calculation shown in Refs.l and 2, Eq.28 is rewritten as 

A - = 8n L j^m\ v ^ (1 " Y T)) ^ (37) 

where is the renormalized velocity of quasiparticles in the superconducting state. See Ref.l 
for the detailed calculation of v^. Y (k; T) is Yosida function: 

Y(k;T) = -Jdel(y^y (38) 



2 

k) ' 



In Eq.38, / (x) = [e x l T + l) 1 is the Fermi distribution function and Ef~ = y s*^ + (A 
where Ajjl = z^A^ is the renormalized superconducting gap. = S a (k,e = 0) is calculated 
by the analytic continuation \e n — ► e + IS from the anomalous self-energy term E a (k). It is 
important that Yosida function is expressed as a universal even function of y = AJj/T: 

g(y)= H e^+^ (l + e^+y^Y 2 dx. (39) 
J —00 ^ ' 



The graph of 1— g (y) for y > is shown in Fig. 2. Note that 1— g (0) = and 1— g (±00) = 1, 
respectively. This indicates that 1—Y (k; T) = at T = T c , since A^ = 0, while 1—Y (k; T) = 1 
at T = 0, and that the increase in y = AJ/T results in the increase in A -2 via Yosida function. 
According to Ref.l, v* (k') = v* (k') at T = T c and v* (k') = j* (k') at T = 0. This means 
that v* (k') is suppressed near T = compared to that near T = T c . From these facts, the 
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i- g{y) 




Fig. 2. The graph of 1 — g (y). For y > 0, 1 — g (y) increases monotonically from zero at y = to 
unity at y = oo. 



inverse square of the penetration depth at T = is given by 

2i f dSfci ^ -* 

X ^\t=0 = ^J fs (2?r) 3 !„• iV^fcV 



(40) 



Here, j£ is the current carried by quasiparticles in the superconducting state and v k is the 
renormalized velocity of quasiparticles in the normal state, respectively. According to Refs.l 
and 2, j'£ coincides with j'£ due to the cancellation of the change of the velocity v k — v k and 
the change of the backflow effect. Therefore we use j k instead of j k hereafter. The velocity of 
quasiparticle v k is expressed as 



d de* 
= z k —(e k + Re^(k,e = 0)) = ~ k 



Jk ~ K dk 

After rather long calculation, 2 j k is expressed as 

fk = < + *kw% (0) , 



dk 



where 



+ 



-T 

k' 



oV 
2vr 



e> Tt f i\ dReGy, (e 1 ) e' dReV. R , {-e') 



2T 



2T 



cV 



(41) 
(42) 

™k (0)=«j?(0) 

^fc'^fc (o), 
(43) 



fc' 



coth ^Im^_ fc , (- £ ') ^^AlZ + tanh — fc -^v lZ ImG g (fi /) 



e> dReVK_ k ,(-e')_ 



de' 



2T 



de' 



(44) 



Zk 



aReSg(fc,e) 
<9e 



£ = 



is the renormalization factor. (fc, e) and the retarded Green's 
function G k (e) is calculated by the analytic continuation \e n — > e + i<5 from S n (k) and G (k), 
respectively. Similarly, V k (e) is derived by the analytic continuation iu;„ — > e + i<5 from the 
normal vertex V (k). 
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Note that = v* (k)\ T=0 is lowered compared to = v* (k)\ T=T . This difference is 
called the vertex correction, which is considered as the backffow effect. According to Ref.22, 
the vertex correction is more effective, when the antiferromagnetic spin fluctuation is strong, 
such as in the underdoped LSCO. The penetration depth \q^ v without including the vertex 
correction is expressed by 



= / t^ttt^m - Y T )) ^ 

Jfs (2vr) \v* k ,\ 

r n = 87F / to * \ v Un v ** ( 46 ) 
T=0 Jfs (27r) \vL\ 



A" 2 



. Since the vertex 

T=0 



Owing to the vertex correction, A~ 2 | T _ is lowered compared to A, 

correction is effective only near T = 0, A~ 2 is suppressed near T = under the strong 
antiferromagnetic fluctuation. In order to calculate AjT^j, at finite temperature < T < T c , we 
have to calculate it in the superconducting state, since Eq.45 contains Yosida function. We 
calculate A^T 2 using Eqs.45, 38 and 41, in which Green's functions and self-energy terms are 
calculated by Eqs. 11-17. Although the vertex correction is not included in Aq 2 , it expresses 
the behavior of A~ 2 near T = T c correctly where the vertex correction is not effective. This 
fact means that we can discuss the development of A~ 2 just below T c using the result of A^,. 

6. The Pressure Dependence of the Parameters 

In this section, we discuss how the parameters in our model depend on the pressure. We 
assume that the dependence of the Coulomb repulsion U on the pressure is weaker than those 
of the hopping terms t and t' . On the other hand, the absolute values of t and t' increase under 
high pressure, since the lattice constant of k-(ET)2Cu(NCS)2 decreases due to compression. 
Owing to the large isothermal compressibility, 23 k-(ET)2Cu(NCS)2 is quite sensitive to the 
pressure. Since the bandwidth W is approximately proportional to t, the large value of t 
results in small ratio of U/W which indicates the weak electron correlation. 

Moreover, we need to take account of the change of t'/t under pressure. In Refs.5 and 6, 
t'/t and the effective mass ratio for a-pocket electron m* a /m e under pressure are measured by 
Shubnikov-de-Haas effect in both K-(d 8 ET) 2 Cu(NCS)2 and k-(/i 8 ET) 2 Cu(NCS) 2 . According 
to these results, the decrease in T c , that in m*/m e , and the increase in t'/t are observed as 
the hydrostatic pressure increases. The decrease in T c is similar to the case of other organic 
superconductors. 24 The decrease in the effective mass ratio is observed also by the infrared 
reflectivity experiment. 25 The stronger pressure dependence of t'/t, T c and m* a /m e in k- 
(d 8 ET) 2 Cu(NCS)2 5 compared to that in k-(/i 8 ET) 2 Cu(NCS) 2 6 is considered to result from 
the different pressure medium used in the two experiments. In fact, the isotope effect for T c 
is small, when the same pressure medium is used. 7 We can consider that T c under pressure 
depends strongly on the pressure medium. 

In K-(d 8 ET) 2 Cu(NCS)2, the ratio t'/t varies from about 0.7 at P = to about 0.77 at 
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P = O.lGPa. The fact that t'/t increases under pressure indicates that t! is more sensitive to 
the pressure than t. On the other hand, m* a /m e varies from 3.5 at P = to 2.5 at P = O.lGPa. 
The change of t'/t results in the change of the Fermi surface. When t'/t is large, the lattice 
structure is near to the triangular lattice, and the Fermi surface is far from the antiferromag- 
netic Brillouin zone boundary as shown in Fig. 8. Therefore the antiferromagnetic fluctuation 
becomes weaker as t'/t increases with the fixed value of U. From these reasons, we can say 
that the electron correlation and the antiferromagnetic fluctuation in k-(ET)2Cu(NCS)2 are 
strong under low pressure compared with those under high pressure. In our model, we as- 
sume that the value of U is independent of the pressure and that the values of t/U and t'/t 
vary with the pressure. From now on, we take U as the unit of energy, since it is assumed 
to be unchanged. In our model, W and the strength of electron correlation mainly depend 
on t/U. On the other hand, the shape of the Fermi surface is changed by t'/t. Note that the 
antiferromagnetic fluctuation becomes weaker as t'/t increases with the constant value of U. 

7. Result and Discussion 

In our numerical calculation, we divide the first Brillouin zone into 128 x 128 meshes in 
the fc-space and take 4096 Matsubara frequencies. In this condition, the obtained results are 
reliable with temperature T down to approximately O.OOlt. Within the FLEX approximation, 
we calculate v^, j kx , Aq 2 | t=0 and A~ 2 | T=Q by using Eqs.41, 42-44, 46 and 40 respectively. 
The results are shown in Figs. 3-7. From now on, we calculate only A~ 2 and refer to it as the 
in-plane penetration depth A~ 2 , since the qualitative behavior of A~ 2 is similar in all directions 
within the xy-plane. We abbreviate the subscripts xx and denote A~ 2 as A -2 hereafter. In the 
calculations of Aq 2 | t=0 and A -2 | T=Q , the temperature at which the calculation is performed 
is important. In order to calculate Eqs.40 and 46 precisely, we need to obtain the values of v kx 
and j kx at T = by extrapolation of v* kx and j kx obtained at finite temperatures. However, 
the extrapolation of j kx is difficult, since the numerical error in the calculation increases in 
the low temperature region. Therefore v kx , j kx , Aq 2 | t=q and A~ 2 | r=0 shown in Figs. 3-5 and 
7 are calculated with the temperature fixed to T = 0.0014^7. Although the calculations are 
not correct, it is enough to investigate the qualitative behavior of the vertex correction. 

Figure 3 shows that j kx is lowered compared to v% x by the vertex correction within the 
FLEX approximation. Note that the ratio of the vertex correction (j£ x — v^ x ) /v\ x is approx- 
imately constant everywhere on the Fermi surface in Fig.3. This is in contrast to that in 
underdoped cuprates calculated in Ref.l. In the cuprates, (j\ x — Vk x )* /v* kx is large at hot 
spots compared to that at cold spots. This implies that in k-(ET)2Cu(NCS)2 the difference 
between hot spots and cold spots (see Fig. 8) on the Fermi surfaces is small in contrast to the 
underdoped cuprates, as long as the temperature near T = is concerned. Figure 4 shows 
z k l , Aq 2 | t=q , and A~ 2 | T=Q for various values of t/U with t'/t = 0.7, where z~ x = z k l is the 
average of z k l over the Fermi surface. The factor z^ 1 is considered to be approximately pro- 
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Fig. 3. The current and the renormalized velocity on the FS for t/U = 0.20, t'/t = 0.7 and the 
temperature T = 0.0014C7 within FLEX. The points 0,P,C1-C4,H1-H4 correspond to those in 
Fig.8. 
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Fig. 4. A~ 2 | T=0 (with VC), A 2 | T=Q (without VC) and z k 1 calculated by FLEX. Low value of t/U 
corresponds to the strong electron correlation and high value to the weak one. the ratio t'/t is 
fixed to 0.7. The vertex correction reduces largely the value of A~ 2 compared to Aq~ 2 , when the 
electron correlation is strong. 



portional to m* a /m e . Figure 5 shows that the effect of the vertex correction within SC-SOPT 
is quite small compared to that in the FLEX shown in Fig. 4. This is similar to the results in 
Ref.2. Figure 6 shows the dependence of T c on t/U with t'/t = 0.7 calculated within FLEX 
and within TOPT. The value of T c increases as t/U decreases both for FLEX and for TOPT. 
For TOPT, small value of t/U is required for superconductivity and T c decreases rapidly as 
t/U increases compared to T c calculated within FLEX. These results are consistent with the 
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calculations in Refs.13, 14 and 15. From Fig.4, we can see that (Aq 2 | t=0 — A _2 | t=q )/ Aq 2 | t=q 
increases as t/U decreases. This means that the suppression of A~ 2 | r=Q by the vertex cor- 
rection is large under the strong electron correlation, even if the shape of the Fermi surface 
is unchanged. Moreover, Aq 2 | t=q itself decreases as t/U decreases. This is because Aq 2 | t=0 
is approximately proportional to zt, or zW, where z = (m*/m)^ 1 is the average of the 
renormalization factor on the Fermi surface. Shortly to say, the strong electron correlation 
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suppresses A _2 | T=0 through the increase of the vertex correction effect and the strong mass 
renormalization which corresponds to the decrease in z. 

Figure 7 shows the t'/t dependence of z^ 1 , Aq 2 | t=q and A~ 2 | T=0 with t/U = 0.200, 
calculated within FLEX. On the other hand, Fig.9 shows the t'/t dependence of T c within 
FLEX and within TOPT. Both in FLEX and in TOPT, T c increases as t'/t decreases and the 
dependence is quite similar to each other. According to Fig.7, (Aq 2 | t=0 — A~ 2 | T=0 )/ Aq 2 | t=0 
increases as t'/t decreases. This means that the effect of the vertex correction in A _2 | T=0 is 
strong under the strong antiferromagnetic fluctuation. Moreover, A~ 2 | T=Q decreases approx- 
imately in proportion to z as t'/t decreases. Then with similar discussion to the dependence 
on t/U, A~ 2 | T=Q is suppressed by the effect of vertex correction and by the decrease in z, 
when t'/t is small. 

Under low pressure, the strong electron correlation and the strong antiferromagnetic fluc- 
tuation exist. In our model, the decrease in t/U and that in t'/t under low pressure indicate the 
strong electron correlation and the strong antiferromagnetic fluctuation, respectively. These 
effects result in the strong mass renormalization and the increase in the vertex correction 
effect. The suppression of A~ 2 | T=Q under low pressure results from these effects. On the other 
hand, the increase of T c under low pressure is also explained by the strong electron correlation 
and the antiferromagnetic fluctuation. 

Let us consider the justification of our calculation by estimating the parameters. By com- 
paring our result for z~ l in Figs. 3-7 with the value of m* a /m e and t'/t measured in Ref.5, we 
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Fig. 8. The Fermi surfaces of ft-(ET) 2 Cu(NCS) 2 with t'/t = 0.7 (dotted line) and 0.75 (solid line), 
where t/U is fixed to 0.200. The shape of the Fermi surfaces is consistent with the estimation from 
the Shubnikov-dc-Haas effect. The crossing points of the FS and the antiferromagnetic Brillouin 
zone boundary shown by the dash dotted line are the hot spots, which are denoted by H1-H4. On 
the other hand, the crossing points of the FS and the gap nodes shown by two diagonal lines are 
the cold spots, which are denoted by C1-C4. Note that the axes are rotated by 7r/4 compared to 
the FS in Refs.5 and 26. 
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Fig. 9. The dependence of T c on t'/t within FLEX and TOPT. The former is calculated with t/U = 
0.200, and the latter with t/U = 0.154. Both for FLEX and TOPT, T c is higher under strong 
antiferromagnetic fluctuations, or small value of t'/t. 



can estimate the parameters in the calculation for K-(<i8ET)2Cu(NCS)2. For convenience, we 
assume z^ 1 to be equivalent to m* a /m e in the estimation. Judging from the values of z" 1 , the 
FLEX approximation is more appropriate for K-((i8ET)2Cu(NCS)2 than SC-SOPT. From this 
reason, we take account of only the FLEX approximation in the discussion of A~ 2 hereafter. 
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Then, we estimate t/U = 0.208 at P = from z' 1 = 3.5 and t'/t = 0.7 5 and t/U = 0.237 
at P = O.lGPa from z' 1 = 2.5 and t'/t = 0.77, by using the result within FLEX. From the 
above estimation, we can compare the calculated T c with the measured one. Since the value 
of t at P = is estimated 6 as 0.12eV ~ 1400K from the polarized infrared reflectance mea- 
surements, 2 U is estimated as approximately 6700K. Using these parameters, we calculate 
T c = 2.08 x 1Q- 3 U = UK at P = 0, and T c = 2.3 x W^U = 1.5K at P = O.lGPa within 
FLEX. The value of T c measured by experiment for K-(dsET)2Cu(NCS)2 is approximately 
10K at P = and 6K at P = O.lGPa ' 28 although it depends strongly on the pressure 
medium. Considering that the calculated T c has delicate dependence on the parameters, we 
can say that the calculated T c is in magnitude consistent with the measured T c and that our 
model and calculation are justified. 

Next, we discuss the development of A~ 2 just below T c in K-(d8ET)2Cu(NCS)2- We cal- 
culate Aq 2 for < T < T c according to Eq.45. As a matter of fact, Aq 2 develops rapidly 
just below T c , when t/U and t'/t are small, as shown in Figs. 10 and 11. Figure 10 shows 
the A 2 -T diagram with t/U = 0.222, t'/t = 0.7 and with t/U = 0.200, t'/t = 0.7. On the 
other hand, Fig.ll shows the diagram with t/U = 0.200, t'/t = 0.7 and with t/U = 0.200, 
t'/t = 0.75. Although the effect of vertex correction is not included in Ag 2 , these diagrams 
express qualitatively correct development of A~ 2 just below T c , since A 2 = A~ 2 at T = T c }' 2 
From Figs. 10 and 11, we can see that Aq 2 has a linear dependence on T near T = 0: 

V 2 - Ao; 2 | T=0 ocT, (47) 

where Aq* | T=0 is obtained from extrapolation of Ag 2 calculated by Eq.45 to T = 0. The 
linear dependence of Aq 2 results from the existence of the line-node in the superconducting 
gap. Ao* 2 | T=0 does not coincide with Ag 2 | T=0 shown in Figs. 4 and 7. This is because the 
results in Figs. 4 and 7 are calculated with T = 0.0014C7 and are not extrapolated to T = 0, 
owing to the difficulties in calculating j kx just near T = 0. Since v kx can be calculated just 
near T = 0, we can extrapolate Eq.46 to T = 0. The difference between A ~ !| 2 | T= g and the 
result of extrapolation of Eq.46 is less than 3%. This indicates that the calculation of Eq.45 
and that of Eq.46 are consistent with each other. 

For the convenience in comparison, we normalize these diagrams by T c and Xq^\ t _ q . The 
results of the normalization for Figs. 10 and 11 are shown by Figs. 12 and 13, respectively. In 
Figs.12 and 13, the development of Ag 2 just below T c with t/U = 0.200 and t'/t = 0.7 is more 
rapid than that with t/U = 0.222 and t'/t = 0.7, or than that with t/U = 0.200 and t'/t = 
0.75. This indicates that the strong electron correlation and the strong antiferromagnetic 
fluctuation are the main origin of the rapid development of A -2 just below T c under low 
pressure in Ref.3. 

The rapid development of Aq 2 just below T c implies that A* k /T develops rapidly. The 
temperature dependence of A femax /T is shown in Fig.14, where AJl max is the maximum of A fe 
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Fig. 10. A 2 -T diagram for (t/U = 0.222, t'/t = 0.7) and for (t/U = 0.200, t'/t = 0.7). T c = 
1.63 x 10 3 C7 in the former, and 2.37 x 10 3 ?7 in the latter. 



on the Fermi surface. In Fig. 14, A* kmax /T develops rapidly with t/U = 0.200 and t'/t = 0.700, 
compared to that with t/U = 0.222 and t'/t = 0.700, or to that with t/U = 0.200 and t'/t = 
0.750. This indicates that the strong antiferromagnetic fluctuation and the strong electron 
correlation cause the rapid development of A£ max /T, which result in the rapid development 
of Aq 2 . Figures 12, 13, and 14 justify the conclusion that the increases in AJj/T is the dominant 
cause of the development of Aq 2 . The pseudogap effect is also considered to be the cause of the 
rapid development of A£ just below T c , although it can not be calculated in the framework of 
this paper. The pseudogap arises from the superconducting fluctuation as discussed by Jujo, et 
al. 29 In fact, the pseudogap effect in rc-(d 8 ET) 2 Cu(NCS) 2 is observed up to T = 45K in STM 
spectroscopy at ambient pressure. 30 ' 31 Although similar measurements under pressure are not 
performed, it is natural to consider that the pseudogap effect is weak or disappears under 
pressure owing to the weak electron correlation and the weak antiferromagnetic fluctuation. If 
the pseudogap exists above T c , then the superconducting gap A£ develops rapidly when the 
superconducting state is realized at T = T c as shown in Ref.32. 

Last, we discuss the uniaxial pressure dependence of T c . We are not able to calculate 
the effect of the interlayer pressure applied parallel to a-axis, since we adopt Q2D band 
structure and ignore the interlayer hopping term. Nevertheless, we can presume that the 
three-dimensionality is strong under the interlayer pressure and it is generally understood 
that T c is suppressed under strong three-dimensionality as calculated in Refs.33 and 34. For 
the 6-axis pressure, the behavior of T c is easily explained. Under the 6-axis pressure, both t 
and t! increase. Since the increase in t is small and that in t' is large under the 6-axis pressure 
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Fig. 12. A 2 -T diagram for (t/f/ = 0.222, t'/t = 0.7) and for (t/U = 0.200, t'/t = 0.7). This diagram 
is normalized by T — T c and Aq 2 = Aq 2 | t _ q , respectively. Aq" 2 shows rapid development under 
the strong electron correlation. (t/U = 0.200) 



compared to those under the hydrostatic pressure, it is sure that both t'/t and t/U increase 
under the 6-axis pressure. Therefore we can say that the 6-axis pressure suppresses T c similarly 
to the hydrostatic pressure from the discussion above. 

Then we discuss the behavior of T c under the c-axis pressure. Since the c-axis is perpendic- 
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0.222 and t'/t = 0.700, and with t/U = 0.200 and t'/t = 0.750. A* kmax /T with t/f/ = 0.200 and 
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ular to the direction of the hopping term t', it is natural to assume that t' is constant and that 
t increases under the c-axis pressure. Moreover, U is assumed to be constant with increasing 
c-axis pressure. These assumptions indicate that t/U increases and t'/t decreases. The former 
corresponds to the suppression of the electron correlations and the latter to the deformation 
of the Fermi surface, which approaches the antiferromagnetic Brillouin zone boundary. These 
two effects are competitive to each other for T c : the former decreases T c and the latter increases 
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it. 

Figures 15 and 16 show the behavior of T c under the c-axis pressure calculated within 
FLEX and TOPT, respectively. For FLEX, T c tends to increase rapidly up to t/U = 0.28 and 
the decrease over t/U = 0.28 is rather slow. The value of t/U = 0.28 corresponds to rather 
higher pressure than lkbar, considering that the increase of t under the uniaxial pressure 
is smaller than that under the hydrostatic pressure. As a conclusion, the value of T c within 
FLEX increases as the c-axis pressure increases, since the deformation of the Fermi surface 
affects strongly compared to the suppression of the electron correlation. This behavior of 
T c is consistent with the results in Refs.8 and 9 below lkbar. This fact implies that the 
antiferromagnetic fluctuation is strong under the weak c-axis pressure. On the other hand, T c 
calculated by TOPT increases slightly up to t/U = 0.16 and then decreases rather rapidly. 
The difference between FLEX and TOPT is mainly due to the t/U dependences of T c shown 
in Fig. 6. T c by TOPT depends rather strongly on t/U than that by FLEX, while the t' jt 
dependences differ little. Although it is quite difficult to estimate the value of t/U for TOPT, 
we estimate t/U at ambient pressure to be from 0.15 to 0.16 judging from the value of T c ~ 11K 
and t ~ 1400K. Therefore, T c calculated within TOPT decreases under the c-axis pressure, 
since the suppression of the electron correlation is dominant over the deformation of the Fermi 
surface. This behavior of T c is consistent with the experimental results 8 ' 9 on T c under the c- 
axis pressure above lkbar. It implies that the antiferromagnetic fluctuation is weak under 
the high c-axis pressure. As a conclusion, we can explain the behavior of T c under the c-axis 
pressure by considering two competitive effects: the suppression of the electron correlation 
and the approach of the Fermi surface to the antiferromagnetic Brillouin zone boundary. In 
the low c-axis pressure region, the latter effect is dominant and FLEX approximation is valid. 
On the other hand, in the high pressure region the former effect is dominant and TOPT is 
valid. 

8. Conclusion 

Under low pressure, the development of A~ 2 just below T c in k-(ET)2Cu(NCS)2 is rapid 
and A -2 near T = is suppressed, when they are compared to those under high pressure. 
These effects, added to the increase of T c , result from the strong electron correlation and the 
strong antiferromagnetic fluctuation under low pressure. This is due to the decrease in t/U 
and that in t'/t. In cuprates, the Fermi surface changes with the change of doping, and the 
antiferromagnetic spin fluctuation is enhanced in underdoped cuprates, while the renormal- 
ization factor z and the bandwidth depend weakly on the doping. 1 ' 2 Then the suppression 
of A -2 near T = in underdoped cuprates are mainly due to the vertex correction. On the 
other hand, the explanation of the pressure dependence of A -2 in k-(ET)2Cu(NCS)2 is rather 
complicated. Under low pressure, the narrow bandwidth W and the strong antiferromagnetic 
fluctuation are observed. The former effect results in the strong electron correlation and the 
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Fig. 15. The c-axis pressure dependence of T c calculated within FLEX. The value of t'/U is fixed to 
0.14 (which has higher T c ) and 0.15(lower). T c tends to increase rapidly up to t/U becomes 0.28. 
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Fig. 16. The c-axis pressure dependence of T c calculated within TOPT. The value of t'/U is fixed 
to 0.100 (which has higher T c ) and 0.108(lower). T c slightly increases up to t/U = 0.16 and then 
decreases rapidly. 



decrease in z, while the latter enhances the vertex correction. In short, the decrease in z and 
the effect of the vertex correction are all the origins of the suppression of A~ 2 at T = 0. 

The rapid development of A~ 2 just below T c under low pressure in k-(ET)2Cu(NCS)2 
is explained by the rapid development of the superconducting gap. The superconducting 
gap develops rapidly just below T c under low pressure, since the strong antiferromagnetic 
fluctuation and the strong electron correlation exist. Moreover, the pseudogap effect is also 
considered to cause the rapid development of the superconducting gap, although it is not 
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studied in this paper. This is similar to the case of cuprates. In cuprates, the strong pseudogap 
effect is observed in the underdoped region where the antiferromagnetic fluctuation is strong, 
and A -2 develops rapidly just below T c compared to that in the overdoped region. 

As a conclusion, we have explained various experimental results on k-(ET)2Cu(NCS)2 un- 
der pressure, e.g. A -2 near T = 0, the development of A -2 just below T c , the superconducting 
transition temperature under the hydrostatic pressure and that under the uniaxial pressure. 
In the explanation, we adopted Hubbard model and assumed that the electron correlation and 
the antiferromagnetic fluctuation are suppressed under pressure. Then the obtained results 
are consistent with the measurements, which indicates that our assumption is appropriate. 

The numerical calculation was performed by sx5 in Yukawa Institute Computer Facility, 
Kyoto University. 
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